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ABSTRACT 


An  analysis  is  carried  out  of  the  effects  of  various  non-ideal  circum¬ 
stances  on  the  two-tone  method  of  measuring  the  two-frequency  correlation 
function  R(f2,  t)  for  microwave  line-of- sight  communication  channels.  These 
circumstances  are  the  presence  of  both  short-term  and  long-term  instabili¬ 
ties  of  the  signal  sources,  the  presence  of  additive  noise,  and  the  use  of 
discontinuous  tones  such  as  those  available  from  a  spinning  satellite. 

It  is  shown  that  oscillator  instability  effects  are  removable  by  assigning 
the  short-term  carrier  oscillator  instabilities  to  the  channel,  using  high 
quality  modulating  oscillators  and  including  phase-lock  tracking  loops,  that 
the  effect  of  additive  noise  can  be  made  negligible  provided  the  signal- to- 
noise  power  ratio  of  the  received  tones  is  greater  than  about  6  db,  and  that 
the  use  of  discontinuous  slowly- switched  tones  is  distinctly  disadvantageous. 

A  method  of  processing  the  received  tones  is  suggested  which  automati¬ 
cally  removes  the  large  constant  delay  due  to  range,  leaving  only  the 
frequency  dependent  part  of  the  delay. 
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The  Two-Frequency  Measurement  of  Line- Of- Sight  Spread  Channels 


I.  THE  TWO-FREQUENCY  METHOD 

In  order  to  transmit  at  a  given  information  rate  with  less  than  a  given 
error  rate  through  a  time-varying  dispersive  channel  it  is  necessary,  in 
general,  that  the  signal  element  be  carefully  chosen  and  various  diversity 
techniques  employed.  This  aspect  of  the  design  of  a  communication  link  is 
made  easier  when  a  general  channel  characterization  is  known  which  is  inde¬ 
pendent  of  the  particular  signal.  One  such  characterization  is  the  two- 
frequency  correlation  function  R(Q,  t)  defined  as 

R(S2,  r)  =  E{H*(f-|,  t-|)  H(f  +  t  +^)}  ,  (1) 

where  H(f,  t)  is  the  time  varying  transfer  function,  E  denotes  ensemble  aver¬ 
age  and  the  star  denotes  complex  conjugate. 

In  general,  R(f2,  r)  is  also  a  function  of  f  and  t,  but  for  line  -  of- sight  at¬ 
mospheric  communication  paths  in  the  microwave  region  it  is  certainly  con¬ 
venient  and  yet  not  too  restrictive  to  assume  that  R(f2,  r)  is  very  weakly 
dependent  upon  f  and  t  compared  with  its  dependence  on  and  r.  This  means 
that  when  R(f2,  t)  is  measured,  it  is  necessary  only  to  specify  the  band  of  the 
microwave  spectrum  in  which  f  lies  and  the  observable  weather  conditions 
that  prevailed  at  the  time  of  the  measurement. 

A  direct  method  of  measuring  R(f2,  t)  is  to  transmit  two  tones. 
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*n(t)  -  exp{ jZiTf^t},  (n  -  1,  2)  over  the  path.  These  are  separated  from  one 
another  at  the  receiver  by  filtering  and  then  time-correlated  with  frequency 
offset  f  —  f^.  The  result  is 

-j2ir(f  -f.)t  ^  jir(f  2  +  f.)T 

<  yl  ^  ”  2>  y2^  +  l)  e  >  =  <H*(f  t  -  j)  H(f2>  t  +  j)  >  e 

jir(f  +  f  )  t 

=  R(f2  ~  fj,  t)  e  ,  (2) 

where  the  further  assumption  of  ergodicity  has  been  made  to  enable  the  time 
average,  denoted  by  <>,  to  be  replaced  by  expected  value.  This  result  de¬ 
rives  directly  from  the  definition  of  H(f,  t),  which  is  that  H(f,  t)  exp{j27Tft}  is 
the  received  signal  when  a  signal  exp{j2Trft}  is  transmitted. 

First,  we  observe  that  the  transmitted  signal  is  not  a  pure  tone  and  so 
we  write 

j27tf  t 

x  (t)  =  a  (t)  e  n  (3) 

n  n 

where  f  is  the  nominal  frequency  and  a  (t)  describes  the  net  effect  of  all 
n  n 

possible  modulating  disturbances  e.  g.  ,  frequency  offset  and  drift,  random 
frequency  jitter,  switching  functions,  if  any,  and  Doppler  effects. 

The  received  signal  y  (t)  is  given  by 


2 


y  (t)  =  l  \  H(f,  t)  X  (t' ; 
n  J  ,]  n 


ej2wf(t-t')dt,df  _ 


(4) 


and  after  mixing  with  b  (t)  exp{j27Tf  t}  to  reduce  the  nominal  frequency  from 


n 


n 


f  to  zero,  the  resulting  signal  z  (t)  is  given  by 
n  n 


z  (t)  =  y  (t)  b  (t)  , 

n  n  n 


(5) 


where  b  (t)  describes  the  net  modulation  of  the  local  oscillator.  (It  is  probably 
more  convenient  in  practice  to  reduce  the  frequency  to  some  common  fre¬ 
quency  different  from  zero,  but  it  is  certainly  possible,  by  quadrature 
processing,  to  choose  a  common  frequency  of  zero  and  it  is  more  convenient 
analytically.  ) 

According  to  Eq.  (2),  R(f^  f  ,  t)  is  measured  by  determining  the  time 

average  of  z  (t  —  r/2)  z  (t  +  t/2).  An  expression  for  this  average  may  be 
n  m 

obtained  by  combining  Eqs.  (3),  (4)  and  (5)  and  then  carrying  out  a  series  of 

transformations.  A  simpler  way  is  to  observe  that  since  x  (t)  is  a  narrow 

n 

band  signal  and  since  the  channel  is  a  microwave  line-of-sight  communication 

path,  then  over  the  band  occupied  by  x  (t),  H(f,  t)  can  be  written  as 

n 

H(f  ,  t)  exp{—  j2iTr  (f  —  f  )/c},  where  r  is  the  effective  range  at  frequency  f 
n  n  n  n  n 

and  c  is  the  speed  of  light.  Then,  from  Eq.  (3)  and  Eq.  (4),  we  obtain 

r  j2rrf  t 

y  (t)  =  H(f  ,  t)  a  (t-— )  e  ”  ,  (6) 

n  n  n  c 
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which,  together  with  Eq.  (5),  yields  the  equation 

<  z  (t  —  ^)  z  (t  +  ?-)  > 
n  Z  m  Z 


=  <HV(f  ,t-^)H(f  ,t  +^)  a  *(t-  —  -^) 
n  Z  m  Z  n  c  Z 


a  (t-—  +h  b  (t-h  b 
m  c  Z  n  Z  m 


(7) 

It  is  quite  reasonable  to  assume  that  the  statistics  of  the  terminals  and 
of  the  channel  are  independent,  so  that  the  time  average  of  the  product  can  be 
written  as  the  product  of  the  time  averages.  In  addition,  from  the  assumption 

T  T 

of  ergodicity  for  the  channel,  <H  (f  ,  t  —  — )  H(f  ,  t  +  — )  >  can  be  written  as 

n  Z  m  Z 

R(f  —  f  ,t).  Thus  Eq.  (7)  reduces  to 
m  n 

<  z  *(t  —  — )  z  (t  +  rr)  >  =  R(f  —  f  ,  t)  <F  (t,  t)  >  (8) 

n  Z  m  Z  m  n  nm 


where 


r  r 

F  (t,  t)  =  a  *(t-  —  -h  a  (t-—  +  hb(t-hb  *(t+7) 
nm  n  cZm  cZnZm  Z 


Now  each  ’’modulating  function,  ”  a  (t)  or  b  (t),  can  be  factored  as 

n  n 


a  (t)  =  a  (t)  A  (t)  , 

n  n  n 


b  (t)  =  J3  (t)  B  (t)  , 

n  n  n 


where  a  ,  B  are  the  ’’random”  parts  of  a  ,  b  and  A  ,  B  the  ’’deterministic” 
n  n  n  n  n  n 


parts.  The  intention  here  is  to  separate  the  short-term  oscillator  instabilities 
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and  other  jitter  effects  from  frequency  offset,  long-term  oscillator  drift  and 


switching  functions.  The  ultimate  specification  of  the  factoring  process  is 
defined  by  the  manner  in  which  the  measuring  equipment  distinguishes  the 
factors. 

We  assume  that  a  (t)  and  S,  (t)  are  sample  functions  from  independent 
n  b 

stationary  ergodic  stochastic  processes.  Then  Eqs.  (8)  and  (9)  can  be 
written  as 


<  z  (t—rr)z  (t  =  R(f  —  f  ,  t)R^  (t)  R  ^‘\t)  <G(t,  t)  > 

n  2  m2  m  n  nm  nm 


(10) 


(a) 

where  R  (r)  is  the  complex  cros  s  -  correlation  function  between  y  (t)  and 


nm 


n 


Y  (t)#  (Y  =  Of,  1 3),  and  G(t,  t)  is  given  by 


G(t,  t)  =  A  *(t-  —  -j)  A  (t-  —  +  j)  B  (t-^)  B  V(t+^) 
cZm  cZn  2  m  2 


n 


(ID 


Thus  the  success  of  the  two-frequency  approach  to  the  measurement  of 

the  channel  characteristics  depends  upon  the  extent  to  which  the  short-term 

(cv )  *(3 ) 

effects,  represented  by  the  factors  R  (t)  and  R  1  (t)  in  (10),  contaminate 

3  nm  nm 

the  result,  and  also  upon  the  errors  made  at  the  receiver  in  estimating  the 
deterministic  features  of  the  transmitted  signals,  represented  by  the  factor 
<  G(t,  t)  >. 

II.  SHORT-TERM  INSTABILITIES 

One  possible  method  of  generating  the  tones  is  to  use  a  single  RF  or 
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carrier  oscillator  modulated  with  various  multiples  of  a  single  reference  os¬ 
cillator.  Thus,  if  the  instantaneous  phases  of  the  two  oscillators  be  repre¬ 
sented  by  cj)  (t)  and  4>^(t),  then  o^(t)  is  given  to  within  an  arbitrary  constant 
factor  by 


a  (t)  =  exp{j<t>  (t)  +  jncj>  (t)}  (12) 

n  c  o 

since  the  output  will  be  stabilized  in  amplitude. 

An  indication  of  the  effect  of  random  variations  in  6  (t)  and  6  (t)  is  ob- 

c  o 

tained  by  assuming  that  the  statistics  of  the  phases  are  independent  and 
Gaussian,  then 


Rla)(T)  =  E{/(t-^  (t+£)} 

~ w  c  m  l 


n 


1  ,  .2  2 

—  —  (m-n)  <r 

R  (t){R  (t)}™  e  ° 

c  o 


(13) 


where  R  (t)  and  R^(t)  are  the  auto-correlation  functions  of  the  output  signals 

2 

from  the  carrier  and  reference  oscillators  respectively,  c is  the  variance 

of  the  reference  oscillator  phase,  and  the  mean  phases  have  been  absorbed  in 

the  arbitrary  constant  factor.  (The  notation  {Ro(t)}mn  has  the  conventional 

interpretation,  viz.  ,  R  (t)  raised  to  the  power  m  times  n.  ) 

A  similar  expression  is  obtained  for  R^  (t)  when  the  same  assumptions 

nm 

are  made  for  S  (t). 

n 

Equation  (13)  demonstrates  the  two  ways  in  which  the  effects  of  oscillator 
instability  is  manifested.  The  first  is  the  reduction  in  the  range  of  t  for 
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which  R(f  —  f  ,  t)  can  be  measured  (high  quality  microwave  sources  have 
m  n 

spectral  widths  typically  of  the  order  of  1  Hz,  so  that  R  (t)  will  be  essentially 

c 

zero  for  |t1  greater  than  a  few  tenths  of  a  second).  The  second  is  the  reduc¬ 
tion  in  the  range  of  f  —  f  for  which  R(f  —  f  ,  t)  can  be  measured,  although 

m  n  m  n 

this  does  not  appear  to  be  serious  for  values  of  f  —  f  of  importance  in  com- 

m  n 

munication  links.  (For  example,  if  a  value  of  f  —  f  of  100  MHz  is  attained 

m  n 

by  multiplying  up  the  signal  from  a  5  MHz  reference  oscillator,  then  the 

2  2 

factor  exp{-(m  -  n)  c r  /2i  is  not  less  than  0.9  provided  c r  is  not  greater 

o  o 

than  about  25  milli- radians,  which  is  well  within  the  capability  of  high  quality 
HF  sources.  ) 

III.  DETERMINISTIC  EFFEC TS 

There  are  two  types  of  deterministic  effects  of  concern  here.  One  in¬ 
volves  the  errors  made  at  the  receiver  in  estimating  the  frequency  and  the 
rate  of  drift  of  frequency  of  the  tones,  and  the  other  is  concerned  with  the 
limitations  imposed  on  the  measurement  by  the  fact  that  each  tone  may  be 
present  for  only  a  part  of  the  time. 

If  the  tone  frequency  and  drift  rate  are  f  and  a  respectively,  and  the 

estimates  of  the  receiver  of  these  quantities  are  f  1  and  a  *  respectively  then 

n  n 

clearly  unless  f  1  —  f  1  =  f  —  f  and  a  ’  —  a  1  =  a  —a  exactly,  the  time 
mnmn  mnmn 

average  of  G(t,  t)  [see  Eq.  (11)]  is  strictly  zero.  In  practice,  of  course,  a 
finite  time  of  averaging  is  used,  so  that  some  error  can  be  tolerated.  The 
obvious  way  to  ensure  that  the  estimates  are  accurate  is  to  employ  phase- 
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lock  tracking,  with  large  time  constant,  at  the  receiver  to  derive  a  "detrended" 
version  of  each  received  tone. 

If  each  "tone"  is  a  slowly  switched  pulse  train  rather  than  being  continu¬ 
ous,  and  the  inter-pulse  period  is  greater  than  the  coherence  time  of  the 
channel  fluctuations,  then  phase-lock  tracking  is  impossible.  In  this  case  a 
possibility  is  to  process  the  signals  by  carrying  out  the  finite  time  average 

using  a  range  of  values  of  f  1  —  f  1  and  a  1  —  a  1  and  then  taking  the  maximum 

m  n  m  n 

of  the  resulting  surface  as  the  required  estimate.  This  two  dimensional  pro¬ 
cessing  must  then  be  repeated  for  each  measuring  point  in  the  f  —  f^,  t 
plane  which  entails  a  large  amount  of  signal  processing.  However,  since  this 
searching  operation  is  effectively  just  another  way  of  phase -lock  tracking, 
then  if  the  one  is  impossible  it  is  unlikely  that  the  other  will  succeed. 

Thus,  for  a  successful  measurement,  it  appears  that  the  switching  period 
of  the  tone  samples  or  pulses  must  be  less  than  the  coherence  time  of  the 
channel  fluctuations.  At  the  even  higher  switching  rate  of  more  than  twice 
the  bandwidth  of  the  channel  fluctuations,  a  single  Fourier  component  of  the 
received  signal  can  be  selected  by  filtering  and  then  processed  as  a  separate 
continuous  tone. 

There  are  three  other  penalties  imposed  by  the  use  of  slowly  switched 
tones,  even  if  the  rate  is  high  enough  for  phase-lock  tracking.  These  are  (a) 
that  the  tone  samples  must  be  tracked  with  a  time  gate,  (in  order  that  the 
maximum  signal- to-noise  ratio  be  obtained),  (b)  that  the  resulting  measure- 
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ment  of  R(f  —  f  ,  t)  is  possible  only  at  intervals  of  t  equal  to  the  switching 
rate,  and  (c),  that  the  averaging  time  must  be  increased  if  the  same  variance 
in  the  measurement  is  to  be  maintained.  Cases  (a)  and  (b)  are  self  evident. 
Case  (c)  is  dealt  with  in  the  next  section. 

IV.  EFFECT  OF  NOISE  AND  FINITE  AVERAGING  TIME 

In  this  section  we  assume  that  the  deterministic  variations  are  not  pre¬ 
sent  except  for  simple  on  -  off  periodic  switching  functions,  that  the  short¬ 
term  equipment  instabilities  are  negligible  and  that  now  additive  white  noise 
is  present.  Thus,  when  a  switched  tone  of  amplitude  a  is  transmitted,  the 
received  signal  is,  after  gating  and  reduction  to  zero  frequency, 

u  (t)  =  s  (t)  {a  H(f  ,  t)  +  w  (t)}  (14) 

n  non  n 

where  s  (t)  is  the  slow-rate  periodic  gating  function  (which  is  synchronized 
n 

with  and  of  shorter  duration  than  the  switching  function  in  the  transmitter,  so 

that  the  latter  does  not  appear)  and  w  (t)  exp{j2nf  t}  is  the  additive  white 

n  n 

noise  associated  with  the  n'th  tone. 

The  measurement  depends  upon  the  finite  time  average 

T_ 

2 

U  (t)  =  -jp  (  u  *(t  —  h  U  (t  +  h  dt  .  (15) 

nm  T  J  n  2  m  2 

T 
“  2 

x 

If  we  focus  attention  on  those  values  of  t  (say  t.)  for  which  s  (t  —  — )  and 

l  n  2 
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*T 

s  (t  +  — )  are  coincident,  then 
m  2 


E{k  AU  (T.)}=a2R(f  -  f  ,  t.)  +  R(w)(t.)  , 

nm  1  o  m  n  1  nm  i 


(16) 


where 


T_ 

Z 

i  r  T ■  T- 

k  =  —  S  (t - s  (t  +  — ^)  dt 

rp  j  n  Z  m  Z 

T 


and 


(w). 


R'  "  (t)  =  E{w  (t  —  rr)  w  (t  +  ^)} 
nm  n  2  m  2 


=  0,  for  n  +  m 


For  our  estimate  of  R(f  *—  f  ,  r.),  therefore,  we  take 

m  n  i 

a  ^{k  U  (t.)  —  RM(t.)},  where  the  noise  contribution  R^W^(t.)  is  differ* 
o  nm  i  nm  l  nm  i 

ent  from  zero  only  when  n  =  m.  It  can  either  be  estimated  directly  by 

measuring  U  (t.)  in  the  absence  of  a  transmitted  signal,  or  be  excluded  by 
nn  i 

basing  the  resulting  plot  of  R(f  —  f  ,  t  )  only  on  those  measurements  of 

m  n 

U  (t)  for  which  n  +  m. 
nm 

Z 

The  variance  cr  (t.)  of  the  estimate  is  computed  in  the  following  manner 


U  (T.)  R  (t.) 

2/  \  t-  nm  1  nm  i  ntr  r  x 

0-  (t.)  =  E-  - - - - - R  f  -  f  ,  T. 

i  12,  2  m  n  l 


(17) 


a  k 
o 
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By  assuming  that  the  channel  statistics  and  the  noise  statistics  are  inde¬ 
pendent  and  that  both  are  multivariate  Gaussian  of  zero  mean,  Eq.  (17)  can 
be  written 


Z 

Z 


o'  (t.)  = 
1 


1 


w 

a  k  T 
o 


n 


a  4|R(0,  f  -  t)!2  +  a  2  R(0,  t’  -  t)  R*(w)  (f  -  t) 
o  i  o  mm 


Z 

2 


+  a  2  R*(0,  f  -  t)  R(w)  (t1  -  t)  +  R(w)  (t1  -  t)  R*(w)  (t'  -  t) 
o  nn  nn  mm 


T ,  T.  T.  T. 

x  s  (t ’  —  — i)  s  (t1  +  —^)  s  (t  —  -r-)  S  (t  +  — -)  dt  dt1 
n  2m  2n2m2 


(18) 


Now  the  four-fold  product  of  the  gating  functions  in  the  integrand  is,  for 

the  special  values  t.,  a  doubly  periodic  array  of  square  columns  standing  on 

the  t,  t'  plane.  If  the  density  of  these  columns  is  high,  that  is,  if  the  period 

of  the  gating  functions  is  short  compared  with  the  width  in  t  of  R(0,  t)  and 

( w ) 

R  (t),  then  the  effect  of  the  gating  on  the  integral  is  merely  to  modify  it  by 

Z 

the  constant  factor  k  .  The  resulting  integral  may  then  be  transformed  to 
give 


CT*  (T 


a  T  T 
o  -  1 


a  4  I R(0,  t)  |2  +  a  2  R(0,  t  )  r'(w)  (t) 
o  I  o  mm 


+  a  2  R*(0,  T)  R  (t)  +R(w)(t)  R¥(w)(t)  }»dt  . 
o  nn  nn  mm 
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Therefore  if  T  is  large  compared  with  the  spread  in  t  of  R(0,  t)  and 

(w )  2 

R  (t),  then  cr  (t.)  is  equal  to  the  sum  of  four  terms  of  the  type 
nn  1  7  r 

^  R  (t)  R?  (t)  dT,  which  is  equal  to  \  S  (f)  S  (f)  df,  where  S  (f)  and  S  (f)  are 

the  spectral  densities  corresponding  to  the  correlation  functions  R^(t)  R^(t). 

The  spectral  density  of  the  noise  can  be  assumed  to  be  N  over  a  band  B, 
the  bandwidth  of  the  receiver  for  each  tone,  and  zero  elsewhere.  The  spec¬ 
tral  density  of  the  channel  variations  may,  for  illustrative  purposes,  be  ap- 

2 

proximated  by  the  Gaussian  shape  A  exp{- 2f  /B  },  so  that  B  is  roughly  the 
2  db  bandwidth  of  the  channel  fluctuations.  Then  since  B  must  be  significantly 
greater  than  B  if  the  tones  are  to  be  processed  without  distortion,  we  obtain 
appr  oximate  ly 


2,  ,  1 
<r  (t.)  ^  — 
l  T 


p2,n  m  2N  R(0,  0)  N  Bn 
R  (0,  0)  +  Q _ +  Q 

<sfir  B  a  2  a  4 


But  a^  R(0,  0)  is  the  received  power  S,  say,  so  that  the  standard  deviation 
normalized  with  respect  to  R(0,  0),  the  maximum  of  R(f2,  t),  can  be  written 
as 


1 


I R(o,  o) 


n/tt 


1  r,  ,  ,  r  V  r  /N\2  Bc 

=  I1  +  2n/*  (s  )  B"+  ^(s  )  B~ 


1/2 


-rr^T  B 


(19) 
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where  S/N  =  S/N^B,  the  signal-to-noise  power  ratio  for  each  received  tone. 

Therefore,  since  B  <  B,  when  the  signal-to-noise  ratio  is  significantly 

greater  than  unity,  the  factor  in  curly  brackets  is  close  to  unity  and  the 

normalized  standard  deviation  of  the  measured  value  of  R(f  —  f  ,  r.)  is 

m  n  i 

roughly  (tt  ^  T  B  )  1  ^ .  Thus,  there  is  little  to  be  gained  by  striving  to  in¬ 

crease  the  signal-to-noise  ratio  above,  say,  £>  db. 

Equation  (19)  has  been  obtained  from  Eq.  (18)  by  using  the  assumption 

that  the  period  of  the  gating  functions  is  short  compared  with  the  width  in 

( w ) 

r  of  R(0,  t)  and  R  (t).  If  we  now  go  to  the  other  extreme  and  assume  that 
nn 

the  period  is  long,  then  only  the  gating  columns  lying  closest  to  the  diagonal 

will  be  effective.  On  the  assumption  that  the  width  of  each  gate  is  small  and 

that  s  (t)  =  s  (t),  we  obtain  in  place  of  Eq.  ( 1  9 ) 
n  m 


<r(r.) 

x 


_1 N 

I R(0»  0)  I  js/TT"  ® 


(19a) 


where  f  is  the  gate  repetition  frequency. 

As  before,  the  optimum  value  of  S  /N  is  of  the  order  of  6  db  but  now  the 

standard  deviation  is  inversely  proportional  to  the  square  root  of  the  number 

(T  f  )  of  samples,  which  is  to  be  expected  since  the  assumption  that  the 

samples  are  independent  is  implicit  in  the  assumption  that  the  switching  rate 

is  slow.  This  assumption  also  means  that  f  <  B  ,  and  therefore,  from  Eqs. 

g  c 

(19)  and  (19a),  that  the  standard  deviation  is  larger  for  the  slowly  switched 
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case. 


The  probable  magnitude  of  B  ,  the  bandwidth  of  the  channel  fluctuations, 

c 

is  of  the  order  of  a  few  Herz.  Therefore,  if  the  normalized  standard  devia¬ 
tion  is  to  be  reduced  to  0.1,  integration  times  of  the  order  of  a  few  tens  of 
seconds  must  be  used.  This  estimate  also  illustrates  the  necessity  of  accu¬ 
rate  frequency  matching  between  the  locally  generated  tones  and  the  received 
tones,  for  the  period  of  the  frequency  match  error  must  be  long  compared 
with  the  integration  time. 

The  average  total  power  output  of  the  transmitter  must  be  divided  be¬ 
tween  the  tones.  It  may  be  advantageous  to  switch  through  the  tones  sequen¬ 
tially  rather  than  radiate  all  tones  continuously.  In  that  case  the  switching 
period  for  each  tone  must  be  less  than  the  Nyquist  interval  of  the  channel 
fluctuations.  This  ensures  that  the  frequency  spread  imposed  by  the  channel 
on  each  Fourier  component  of  the  switched  tones  does  not  overlap  the  spread 
imposed  on  the  neighboring  tones.  Then  the  spread  central  Fourier  com¬ 
ponent  can  be  selected  without  distortion  by  filtering. 

The  question  arises  of  whether  an  improved  signal-to-noise  ratio  can  be 
attained  by  making  the  radiated  tone  samples  very  short  and  of  high  intensity. 
It  is  easy  to  show  that  if  P  is  the  average  power  of  each  tone  and  the  "tone11 
is  radiated  for  a  fraction  1  / m  of  the  time,  then  after  the  central  Fourier 
component  has  been  selected  by  filtering,  the  resulting  continuous  signal  has 
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an  average  power  given  by  P/m.  (The  peak  power  is  mP  but  the  power  in 

2 

the  central  component  is  1/m  of  this.  )  Therefore,  if  N  is  the  additive  white 

o 

noise  power  per  unit  bandwidth,  the  signal- to-noise  power  ratio  is  P/m  N  B 
if  no  pulse  gating  is  used  at  the  receiver  and  is  P/N  B  if  gating  is  used.  Here 
B  is  the  bandwidth  of  the  receiver  for  each  tone. 

The  answer  to  the  question,  therefore,  is  that  no  improvement  is  possible 
by  using  very  short  high  intensity  pulses,  since  the  best  obtainable  signal-to- 
noise  power  ratio  P/N  B  is  independent  of  m. 

V.  PHASE  CONTROL 

Measurement  of  the  two  frequency  correlation  function  involves  the 
measurement  of  both  amplitude  and  phase.  Thus,  it  is  necessary  to  keep 
close  control  of  the  phase  relationships  between  the  tones.  Clearly,  if  each 
tone  is  separately  tracked  by  a  phase-lock  loop,  the  differences  in  phase  be¬ 
tween  the  tones  are  eliminated  and  no  meaningful  phase  measurement  can  be 
made. 

However,  if  a  multiple  tone  transmission  is  used  in  which  the  phase  inter¬ 
relationships  between  the  tones  is  fixed,  it  seems  that  it  is  possible  to  use 
phase-lock  tracking  and  still  retain  the  required  phase  information.  For 
example,  the  transmitted  tones  may  be  generated  by  using  a  single  carrier 
oscillator  mixed  with  various  integral  multiples  of  the  output  from  a  single 
reference  oscillator.  At  the  receiver,  the  necessary  local  oscillator  tones 
would  be  generated  using  an  identical  mixing  procedure  to  that  used  for  gen- 
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erating  the  tones  in  the  transmitter.  The  two  oscillators  in  the  receiver, 
however,  would  be  voltage  controlled  by  phase -lock  tracking  just  two  of  the 
transmitted  tones  with  the  two  corresponding  locally  generated  tones.  Then 
if  the  channel  were  strictly  linear  phase,  the  remaining  tones  would  also  be 
effectively  phase  locked,  but  any  deviation  from  a  linear  phase  law  would 
manifest  itself  as  phase  errors  between  the  remaining  unlocked  tones  and 
their  corresponding  locally  generated  tones. 

This  method  has  the  interesting  property  that  range  effects  are  automati¬ 
cally  taken  into  account  so  that  there  is  no  need  to  estimate  the  expected 
(large)  number  of  phase  rotations  between  one  tone  and  the  next  in  order  to 
obtain  the  particular  phase  interval  in  which  to  place  the  measured  phase.  In 
practical  terms  the  method  has  the  effect  of  removing  the  constant  time  delay 
due  to  range,  which  is  of  no  interest  in  any  case,  and  leaving  only  the  fre¬ 
quency  dependent  part  of  the  delay. 

The  method  is  examined  in  detail  in  the  Appendix  where  it  is  shown  that 
the  two  tones  chosen  for  tracking  must  be  the  two  with  the  smallest  frequency 
difference  between  them.  If  this  is  not  done,  ambiguity  can  arise  in  the  defi¬ 
nition  of  the  phase  of  some  of  the  remaining  tones.  It  is  also  shown  in  the 
Appendix  that  the  full  benefits  of  the  method  cannot  be  attained  when  just  two 
tones  are  transmitted  during  each  time  averaging  interval,  even  if  the  two 
tones  are  derived  from  a  single  oscillator. 

The  fact  that  this  method  requires  the  simultaneous  transmission  of  at 
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least  three  of  the  full  set  of  tones  required  for  the  complete  measurement  is 


not  a  disadvantage  of  the  method  because  the  variability  of  weather  conditions 
makes  it  necessary  to  carry  out  the  complete  measurement  in  as  short  a  time 
as  possible.  Thus  multiple-tone  transmission  is  required  in  any  case. 

It  is  necessary,  of  course,  to  compensate  for  any  difference  in  phase 
rotation  occurring  between  the  various  tone  channels. 

VI.  CONCLUSION 

In  view  of  the  assumption  of  stationarity,  which  is  basic  to  the  charac¬ 
terization  of  the  channel  by  the  two-frequency  correlation  function,  it  is 
important  that  all  the  measurements  necessary  to  plot  the  function  be  made 
simultaneously  and  in  as  short  a  time  as  possible.  The  two-frequency 
method  of  measurement  is  well  suited  to  these  requirements  because  a  large 
number  of  tones  can  be  transmitted  simultaneously  and  compared  in  pairs  to 
cover  the  whole  range  of  frequency  differences  of  interest.  Also,  since  the 
receiver  bandwidth  can  be  made  as  small  as  the  net  uncertainty  in  frequency 
of  the  received  tones,  the  method  can  discriminate  well  against  additive 
noise  and  therefore  make  do  with  a  short  averaging  time  to  achieve  the  de¬ 
sired  variance  in  the  measurements. 

The  facts  that  the  transmitter  may  be  moving  with  respect  to  the  receiver 
and  that  the  frequencies  of  the  transmitted  tones  are  not  known  exactly  at  the 
receiver  make  it  necessary  to  carry  out  some  frequency  estimating  procedure 
at  the  receiver.  The  only  practicable  way  of  dealing  with  this  problem  appears 
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to  be  phase-lock  tracking. 

If  phase-lock  tracking  is  used  on  every  received  tone  phase  information 
is  lost.  However,  it  has  been  shown  by  using  a  suitable  method  of  tone  gen¬ 
eration  and  a  selective  use  of  phase-lock  tracking  all  useful  phase  information 
can  be  measured  directly  since  the  pure-delay  linear  phase  is  cancelled  out. 

An  objection  to  the  choice  of  the  two-frequency  method  of  measurement 
might  be  that  oscillator  instabilities  (and  especially  those  of  the  carrier  os¬ 
cillator)  preclude  any  accurate  measurement  for  values  of  r  other  than  t  =  0. 
However,  it  is  not  possible  to  use  the  channel  without  a  carrier  oscillator  and 
so  it  is  quite  in  order  to  define  the  "channel 11  to  be  the  combination  of  modu¬ 
lation,  transmission  and  demodulation  operations  taken  together.  Then,  since 
the  time  correlation  of  the  channel  is  as  least  as  bad  as  that  of  the  carrier 
oscillator,  it  should  be  easy  to  attain  sufficient  short-term  stability  in  the 
other  lower  frequency  oscillators  involved  in  the  measurement. 

Finally,  the  further  problems  associated  with  the  use  of  slowly  switched 
tones  have  been  examined  in  some  detail.  This  is  because  it  is  necessary  to 
have  either  a  full- stabilized  vehicle  or  an  elaborate  antenna  structure  for  a 
satellite  source  to  appear  non-fluctuating  at  the  ground.  If  the  measurement 
can  be  made  when  the  source  is  single  narrow-beam  antenna  mounted  on  a 
spin- stabilized  satellite  the  experiment  is  much  easier  to  implement. 

The  results  show  that  the  effects  of  using  such  "slowly- switched"  signals 

are: 
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(1)  Phase-lock  tracking  is  no  longer  possible,  which  means  that 
extensive  data  processing  is  necessary  at  each  measured  point 
to  search  for  the  appropriate  frequency  drift  and  offset, 

(2)  Measurements  can  be  made  only  at  values  of  t  equal  to  integral 
multiples  of  the  switching  period, 

(3)  Time-gate  tracking  is  necessary  in  order  to  maintain  the  maxi¬ 
mum  signal- to-noise  ratio, 

(4)  The  averaging  time  necessary  to  achieve  a  given  variance  in 
the  measurements  is  increased. 

Thus  the  use  of  continuous  tones  (or  tones  sampled  at  intervals  no  larger  than 
the  Nyquist  interval)  is  much  to  be  preferred. 
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APPENDIX 


This  appendix  demonstrates  the  possibility  of  using  phase -lock  tracking 
in  such  a  way  that  essential  phase  information  is  retained  and  the  constant 
range  delay  is  automatically  cancelled. 

We  consider  a  set  of  transmitted  tones  generated  by  mixing  multiples  of 
a  reference  oscillator  signal  with  a  carrier  oscillator.  Then  the  phase  of 
each  tone  may  be  written  as 

Arg{x  (t)}  =  cp  (t)  +  ncp  (t)  (A-i) 

n  c  o 

where,  since  we  are  concerned  here  solely  with  the  slowly  varying  or 

"deterministic"  variation,  cp  (t)  and  c p  (t)  are  the  smoothed  instantaneous 

c  o 

phases  of  the  carrier  oscillator  and  the  reference  oscillator,  respectively. 

By  the  time  the  tone  reaches  the  receiver  it  has  received  a  phase  shift 
which  is  approximately  proportional  to  the  frequency,  and  so  the  phase  of  the 
received  tone  can  be  expressed  as 

Arg{y  (t)}  =  cp  (t)  +  ncp  (t)  —  [cp ' (t)  +  ncp ' (t) ]  p  ,  (A-2) 

o  c  o  c  on 

where  the  primes  denote  differentiation  with  respect  to  t  (so  that  <p^(t)  + 
ncp  1  (t)  is  2tt  times  the  smoothed  frequency)  and  is  a  factor  weakly  dependent 
upon  both  the  frequency  and  the  time.  Thus  may  be  regarded  as  the  effec¬ 
tive  frequency  -  dependent  instantaneous  range  divided  by  the  speed  of  light. 

The  suffix  n  denotes  that  p  is  evaluated  at  the  frequency  cp  1  (t)  +  ncp  ’(t). 

n  c  o 

At  the  receiver  a  similar  set  of  tones  is  generated  to  provide  the  local 
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oscillator  signals.  The  phase  of  these  tones  may  be  written  as 

Arg{b  (t)}  =  ip  (t)  +  nip  (t)  ,  (A-3) 

n  c  o 

where  ijj  (t)  and  (t)  are  the  phases  of  the  receiver  oscillators  corresponding 

to  the  carrier  and  reference  oscillators  in  the  transmitter,  respectively. 

The  receiver  oscillators  are  voltage  controlled  by  means  of  phase-lock 

loops  which  lock  the  phases  of  two  of  the  locally  generated  tones  with  the 

smoothed  phases  of  their  corresponding  received  tones.  We  may  choose, 

without  any  loss  of  generality,  to  lock  b  (t)  and  b  (t)  with  y  (t)  and  y  (t). 

°  g  °  g 

That  is,  we  adjust  i|j  (t)  and  ^  (t)  to  satisfy  the  equations 

Arg{b  (t)}  =  Ar g{y  (t) }  , 

o  o 

Arg{b  (t)}  =  Arg{y  (t)}  . 

g  g 

Written  in  full  these  equations  are,  from  (A- 2)  and  (A-3), 
ip  (t)  =  <p  (t)  —  <p ' (t)  p 

C  C  C  O 

(t)  +  gip  (t)  =  <p  (t)  —  <p  '(t)  p  +  g{(p  (t)  -  (p  '(t)  p  }  +  k2tr  . 
c  o  c  c  g  o  O  g 

The  k2iT  term  is  necessary,  where  k  is  any  integer,  because  any  integral 
multiple  of  2tt / g  can  be  added  to  vp  (t )  without  altering  the  conditions  for  phase 
lock. 

By  solving  these  equations  for  ijj  (t)  and  ^o(t)  we  find  that  the  phase  of 
b^(t)  is  given  by 
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gp  +  np  —  np 

Arg{b  (t)}  =  (p  (t)  +  rup  (t)  +  —  2tt  —  ip  '(t)  — - - - —  —  rup ' (t)  p 

neogc  g  og 

(A-4) 

Therefore  when  the  received  tones  are  reduced  to  a  nominal  frequency  of 
zero  by  mixing  with  their  corresponding  locally  generated  tones  given  by  Eq. 
(A-4),  we  find  for  the  phase  of  the  resulting  signals 


Arg{z  (t)}  =  Arg{y  (t)}— Arg{b  (t)} 
n  n  n 


gp  +  zip  —  np 

=  —  {<P ' (t)  +  rup  *  (t) }  p  +  cp  ' (t) - - - - —  +  rup ' (t)  p 

c  o  ne  g  o  g 


(A-5) 


From  Eq.  (A-5)  we  see  that  Arg{z  (t)}  =  0  and  Arg{z  (t)}  =  —  1c2t t,  so 

^  § 

that  the  measured  phase  will  be  zero  for  n  =  0  and  n  =  g,  as  it  should  be 
since  these  values  of  n  are  the  reference  points  of  the  phase  control.  Further, 
if  the  medium  is  non-dispersive  then  p  will  be  independent  of  n  and  in  this 
case  Arg{z  (t ) }  will  be  equal  to  —  rikZir/g .  Thus  if  g  is  chosen  to  be  unity, 
Arg(zo(t)}  will  be  measured  as  zero  in  the  non-dispersive  case  for  all  n. 

Therefore,  we  choose  x  (t)  and  x.  (t)  as  our  reference  tones,  and  Eq.  (A-5) 

o  1 

reduces  to 


Arg{z  (t)}  =  —  {(p ' (t)  +  rup  '(t)}  p  +  <p  '  ( t)  (p  +  np 
n  c  o  n  c  o  1 


np  )  + 
o 


n<p  1  (t)  p 

o  1 


(A-6) 
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That  this  expression  does  in  fact  equal  the  phase  of  z^(t)  relative  to  the 
linear  phase  law  defined  by  the  reference  points  is  easily  tested.  If  we  re¬ 
present  the  phase  difference  over  the  channel  for  the  n'th  tone  as  0  then  0 

n  n 

is  given  by  0^  =  —  {</?  ’(t)  +  mp^(t)}  p  .  The  straight  line  phase  0^  based  on 

0  and  G  .  as  reference  points  is  given  by  0=0+  n(0  .  -  0  ),  so  that  the 
ol  nolo 

relative  phase  difference  AO  between  the  actual  phase  0  and  the  linear  phase 

n 

0  is  given  by  AO  =0  —  0  —  n(0  .  —  0  ).  When  the  full  expressions  for  the 

n  nolo 

0  are  substituted  in  this  equation  it  becomes  identical  to  Eq.  (A-6). 
n 

It  is  easy  to  show  that  if  a  pulse  of  bandwidth  B  is  transmitted  through  a 
dispersive  channel  having  a  quadrature  component  in  the  phase  law,  then  the 
onset  of  severe  distortion  occurs  when  the  deviation  of  the  phase  law  from  a 
linear  law  over  this  bandwidth  is  of  the  order  of  one  radian.  Since  we  are 
interested  in  the  measurement  of  communi cation  channels,  therefore,  no 
ambiguity  in  the  measurement  of  the  phase  should  arise. 

If  the  channel  is  strictly  stationary  it  is  possible  in  principle  to  measure 
the  two-frequency  correlation  function  by  transmitting  a  pair  of  tones  which 
are  changed  in  frequency  after  each  averaging  interval.  By  this  means  a 
better  signal-to-noise  ratio  is  obtained  since  the  total  power  is  divided  be¬ 
tween  only  two  tones. 

Clearly,  a  phase-lock  system  of  the  type  described  above  will  not  work 
in  this  case  because  once  both  tones  have  been  used  for  locking,  all  deter¬ 
ministic  phase  information  is  destroyed.  However,  if  only  one  oscillator  is 
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used  to  generate  both  tones  there  is  just  one  oscillator  to  be  locked  at  the  re¬ 
ceiver  and  the  possibility  arises  of  achieving  the  same  performance  as  in  the 
previous  case. 

The  phase  of  the  received  tones  is  n{c p  (t)  —  (p  1  (t)  p  },  the  phase  of  the 

o  on 

locally  generated  tones  isnij^t)  and  so  the  equation  to  be  satisfied  for  phase 
locking  on  one  of  the  tones  is 

gi}jo(t)  =  g{<PQ(t))  —  <p°  (t)  p  }  +  k2ir  . 

This  equation  gives  i|i  (t)  directly  and  so  the  result  of  mixing  the  other  re¬ 
ceived  tone  with  its  corresponding  locally  generated  tone  is  the  phase 
— *  n(p  1  (t)  (p  — •  — *  nk277/g.  Since  it  is  not  possible  in  this  case  to  choose  g 

equal  to  unity  (in  fact,  g  must  be  a  large  integer)  there  remains  the  irreduc¬ 
ible  phase  uncertainty  nkZir/g.  Thus,  we  conclude  that  when  the  tones  are 
transmitted  two  at  a  time  the  phase  lock  technique  cannot  be  used  to  full  ad¬ 
vantage. 
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